Yukawa Institute Kyoto 



YITP-12-76 



Confluence of apparent singularities in multi-indexed 
orthogonal polynomials: the Jacobi case 



C.-L. Ho ab , R. Sasaki 6 and K. Takemura c 

a Department of Physics, Tamkang University, Tamsui 25137, Taiwan^ 
§"! b Yukawa Institute for Theoretical Physics, Kyoto University, 



C/3 



Kyoto 606-8502, Japan 
c Department of Mathematics, Faculty of Science and Technology, 
Chuo University, 1-13-27 Kasuga, Bunkyo-ku Tokyo 112-8551, Japan 



Abstract 

The multi-indexed Jacobi polynomials are the main part of the eigenfunctions of 
exactly solvable quantum mechanical systems obtained by certain deformations of the 
Poschl- Teller potential (Odake-Sasaki). By fine-tuning the parameter (s) of the Poschl- 
Teller potential, we obtain several families of explicit and global solutions of certain 
second order Fuchsian differential equations with an apparent singularity of character- 
istic exponent —2 and —1. They form orthogonal polynomials over x G (—1,1) with 
weight functions of the form (1 — x) a (l + x)^ /{(ax + b) 4 q(x) 2 }, in which q(x) is a 
polynomial in x. 

1 Introduction 

The theory of Fuchsian differential equations is essentially a local theory [TJ [2] . Exactly solv- 
able quantum mechanics, on the other hand, has recently provided infinitely many examples 
of global solutions, called the multi-indexed Jacobi polynomials j3], of second order Fuchsian 
equations with 3 + I (I > 1) regular singularities. The first three regular singularities corre- 
spond to those for the Jacobi polynomials and the extra £ singularities are all apparent and 
their characteristic exponents are generically all —1. At this point a naive and deep question 
is Hf it is possible to construct global, i.e. rational solutions having higher (p = —2, —3, . . .) 
characteristic exponents at apparent singularities?' In a sense, this could be regarded as 
a problem of the confluence of apparent singularities. In a previous paper [I], a related 
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question was asked in connection with the multi-indexed Laguerre polynomials [3j [5] and a 
partial positive answer was obtained. Namely, several examples of second order differential 
equations with one irregular singularity at oo and two or more regular singularities at finite 
locations were constructed together with solutions having a characteristic exponent —2 at 
one of the regular singularities, through fine-tuning of the parameter of the multi-indexed 
Laguerre polynomials. 

In the present paper, we provide a partial positive answer to the above question within 
the framework of the multi-indexed Jacobi polynomials. As recapitulated in section two, the 
multi-indexed Jacobi polynomials are obtained as the main part of the eigenfunctions of a 
quantum system with the Poschl- Teller potential [HI [7], deformed by Darboux-Crum transfor- 
mations [EJEJCEO] hi terms of the polynomial solutions called the virtual state solutions. The 
"confluence of apparent singularities" is achieved through fine-tuning of the two parameters 
(g, h) of the Jacobi polynomial to a curve g = g(h) or h = h(g) in such a way that the seed 
solutions have a cubic zero. In other words, these curves are the roots of the discriminants of 
the polynomial part of the candidate seed solutions, see (I2.23p . Here we show several exam- 
ples of such Fuchsian equations (with a p = —2) together with their global solutions. We also 
restrict the parameter ranges such that the extra singularities do not appear in the physical 
region. Therefore the fine-tuned multi-indexed Jacobi polynomials form orthogonal polyno- 
mials over (—1, 1) with the weight function of the form (1 — x) a (l + x) l3 /{(ax + b) 4 q(x) 2 }, in 
which q(x) is a polynomial in x. The denominator does not vanish in (—1, 1). Most of the 
explicit examples in section four have four regular singularities at 0, 1, t and oo, for which 
the Heun's equation is the canonical form: 

d 2 y ( 7 5 e \ dy a(3z — q 



dz 2 ' \z z — 1 ~^ z — t)dz^~ z(z — i)(z — t)^ ' ^ ^ 

with the condition r y + 5 + e = a + /3 + l (Fuchs' relation). As is well known (0, 1 — 7), 
(0, 1 — S) and (0, 1 — e) are the exponents at 0, 1 and t, respectively and (a, (3) are those at 
00. The parameter q is called the accessory parameter. Heun's differential equation appears 
in several systems in physics [TTI IT2] , in particular, its special cases appear as the differential 
equations for the exceptional Jacobi polynomials of type Xi, and some of their properties 
were studied in [13J. By the change of the independent variable 77 = 1 — 2z, the equations 
satisfied by the polynomials in all the examples of section four, (14.111) . (I4.12p . flA.4j) . flA.8j) . 
(1A.14I) . (1A.20j) can be cast into the above canonical form (II. ip . 
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This paper is organised as follows. In section two, the essence of the multi-indexed 
Jacobi polynomials is recapitulated. Starting with the general Schrodinger equation with 
Darboux-Crum transformations in section I2.1[ the details of the Poschl- Teller potential and 
the eigenpolynomials as well as the virtual state solutions are introduced in section 12.21 
and the explicit forms of the multi- indexed Jacobi polynomials are explained in section 12.31 
The method to realise the higher degrees of apparent singularities is explained in section 
|3j Section @] is the main body of the paper. In section 14.11 eight distinct seed functions 
are chosen among the Wronskians of two and three virtual state solutions of lower degrees. 
Their parameters are fine-tuned so that a cubic zero is realised. They are divided into four 
groups, I, II, III and IV. The Groups I and II (III and IV) have the same form of potentials 
( 14. 2p . (respectively (I4.4p ) and thus the same form of differential equations (14. 9p and (14. lip 
(respectively (I4.10p and (I4.12p ) when expressed in the shifted parameters g and h (I2.28p . 
Within the same group, the potential and the differential equations of one member's can 
be obtained from the other's by a simple shift of the parameter, e.g. h — > h — 2, etc. In 
section 14.21 the aspects of Heun's equation are discussed in some detail. The limit to the 
multi-indexed Laguerre systems is discussed in section 14.31 Several other interesting cases 
are commented on in section [5j The final section is for a summary and comments. Appendix 
lAl provides the explicit forms of the fine-tuned multi-indexed polynomials, the differential 
equations for the rational as well as the polynomial solutions and the orthonormality relations 
for the Groups I, II, III and IV, together with the Laguerre limits of IV. 

2 Multi- indexed Jacobi polynomials 

Here we recapitulate the essence of the multi- indexed Jacobi polynomials as derived in [3]. 
The main ingredients are (i) the Jacobi polynomials as the main part of the eigenf unctions 
of the exactly solvable quantum system, the Poschl- Teller potential, (ii) the virtual state 
solutions, (iii) the Darboux-Crum transformations. 

2.1 Darboux-Crum transformations 

We start with the Darboux-Crum transformation, which applies to any one-dimensional 
Schrodinger equation. Let 

n = ~ + u(x), (2.i) 
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be a Hamiltonian and U(x) be a potential denned in a complex domain. To define the 
Darboux-Crum transformation, we pick up distinct seed solutions {ipj(x),Sj} (j — 1, . . . , M) 
of the original Schrodinger equation. Namely, 

H<p j (x) = S j <p j (x), SjEC, j = 1,...,M, (2.2) 

and the functions <fj{x) need not be square- integrable. Let 

U {M) 4>{x) = g^x), (2.3) 

ax 2 dx 2 
be the M-th deformed Schrodinger equation, where W[/i, . . . , /jv](a;) is a Wronskian 

W[f 1 ,...J N ](x) d = i Det(^>P-) . (2.5) 



l<j,k<N 



Then we have 

Proposition 2.1 [91 [10] Let ip(x) be a solution of the original Schrodinger equation 

Hil>(x) = £ip(x). (2.6) 

Then the function 

(M){x)d jf W[^,...^ M Mx) 

K } W[ip 1 ,...,i PM ]{x) K ' 

satisfies the M-th deformed Schrodinger equation (12. 3p with the same energy: 

H {M) ^ M \x) = £^ M \x). (2.8) 

In the M-th deformed potential (12.41) . the zeros of the Wronskian W[</?i, . . . , (Pm](x) are 
the extra singularities. 

2.2 Poschl-Teller potential 

The Jacobi polynomials constitute the main part of the eigenfunctions of a quantum system 
with the Poschl-Teller (PT) potential: 

H = -f- 2 +U(x;g,h), U(x;g, h) = 9 -^± + - (g + h)\ (2.9) 

dx 2 sin x cos 2 x 



g>0, h>0, xg(0,|). (2.10) 
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It has regular singularities at x = and x = n/2 modulo tt, and the characteristic exponents 
at x = (resp. x = vr/2) are p — g, 1 — g (resp. p — h, 1 — h). The Hamiltonian is self-adjoint 
and the eigenfunctions and the eigenenergies are given by 

(f> n {x;g,h) = (smx) 9 {cosx) h Pi 9 - l/2 ' h - 1/2 \r](x)), r]{x) = cos(2x), (2.11) 
S n (g,h) = 4n{n + g + h), n = 0,1,2,..., (2.12) 

where P^'^ (if) is the Jacobi polynomial in the variable r] defined by 

p[c,p) {r) \ _ ( a + l)n (-n) k (n + a + (3 + l) k f l-r] \ k 

Pn {v) -^^l< WTTjl {—) ■ (2 ' 13) 

The eigenpolynomial solutions form a complete set of orthogonal polynomials over (0,7r/2) 
under a suitable restriction of the parameters g, h: 

»n/2 



-*n) 0m) 



<fi n (x; g, h)(f) m (x; g, h)dx = h n (g, h)6 nm , (2.14) 

2nl{2n + g + njl (n + g + h) 

Another ingredient is the virtual state solutions, which are three types of seed polynomial 
solutions of the Schrodinger equation (12. 9p indexed by v 6 Z >0 , which have all real energies: 

ft y (x;g,h) = (sinx) 9 (cosx) 1 - /l P v (9 - 1/2 ' 1/2 - fe) (r/(a;)), 

£l(g, h) = -4(g + v + l/2)(h - v - 1/2); (2.16) 
$(x; g, h) = (sinx) 1 - 9 (cosx) h P v (1/2 - 9 ' /l - 1/2) (r7(x)), 

£?(g, h) = -A(g - v - 1/2) (/i + v + 1/2); (2.17) 
^\x;g,h) = (smxy-^cosxy^P^-^-'^i^x)), 

£l u (g, h) = -4(v + + /i - 1 - v) = £- {v+1) (g, h). (2.18) 

Note that the seed solutions of type I and II are used to construct generic multi-indexed 
Jacobi polynomials in [3]. The third type III is now employed to produce a wider class of 
polynomials, to be called tentatively non-generic multi-indexed Jacobi polynomials (for the 
non-generic multi-indexed Laguerre polynomials, see jl]). Because of the type III virtual 
solutions, the absence of singularities in rj G (—1,1) is no longer guaranteed. One salient 
feature of the non-generic multi-indexed polynomials is the existence of the additional modes 
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at n < 0. When 0* n is used as a seed solution, an additional mode at n = — (v+ 1) is created. 
The degree of the polynomial 4>\(x) (resp. <fry(x), is equal to v under the condition 

g-h + k^0 (resp. g - h - k ^ 0, g + h + 1 - k ^ 0) for Jfe = v + 1, v + 2, . . . , 2v. (2.19) 

In the rest of this section, we assume that the parameters are general and satisfy these 
conditions, unless otherwise stated. In section |3] and after, fine tuning of the parameters will 
be used extensively and the above conditions will not be assumed any longer. 

2.3 Explicit forms of the mult i- indexed Jacobi polynomials 

We define a wider class of multi-indexed Jacobi polynomials by applying Darboux-Crum 
transformation with three types of the seed solutions to the quantum system with the PT 
potential in (T279T) . Let . . . , ^ (ti, . . .t/v G {I, II, III}, Vi, . . . Vyy £ Z >0 ) be distinct seed 
solutions. Let us denote the number of type I seed solutions in the set {<p^\, ■ ■ ■ ,0*^} by 
M > 0, that of type II by iV > and that of type III by L > 0. Then Af = M + N + L is the 
total number of the used seed solutions. It follows from Proposition 12.11 that the function 

eKx) de f w[e...,e, 0w] ^ (2 ^ 2o) 

W[$i,...,0$](ar) 
is an eigenfunction of the deformed Hamiltonian 

wW = _^ + ^ ;s , ft) _ 2 «!^W (2 . 21) 

with the same eigenvalue S n = 4n(n + g + h), provided that the deformed potential is 
non-singular. Let us express the numerator and the denominator of <p^\x) as 

W[<, . . . , 0n ]( x ) <^ f (i _ r? )A (1 + rjfQ n ( v ), Q n (±l) ^ 0, (2.22) 
W[<, • • • , $»){x) d = f (1 - + vfV(v), V{±1) ± 0, (2.23) 

in which Q n (v) an d are polynomials in rj. 
Proposition 2.2 T7ie exponents are given by 

A^i(J + l)(2s + J), S d = / i(ir+l)(2/ i + ir), (2.24) 

A' = f ^J(2g + J-l), E' d = f ±K(2h + K- 1), (2.25) 
m& J = M-N-L, K = -M + N-L. (2.26) 
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The deformed potential is given by 

uW(x) *? U-9 ~ 1) 2h(h-l) _ - 2 

W 1 - 7/ 1 + 7/ V ' 



+ 



2 V'(r ] Y-V(r ] )V"(ri) Mrj) 
V) ^fo)J' 



(2.27) 



g = g + M + N-L, h = h — M + N — L, g + h + 2L = g + h } (2.28) 
and iTie eigenfunctions have a simple form 

<t>W(x) = (1 - r,) 5/2 (l + l) i/2 -|^. deg(fi n ) = £ + n + L, deg(P) = C, (2.29) 

( s "g« yj _ M(M-l) _ Af(Af-l) _ L(L-1) + MJV ^ Q p 30) 

2 2 2 

For general values of g and h, the two polynomials Q n {v) an d ^ > ( 7 ?) are relatively prime and 
the polynomials {Q n (v) = 2n(n; g, h)} are called the multi-indexed Jacobi polynomials. 

With a slight change of notation, these formulas are generalisation of those in [3] with the 
inclusion of the type III seed solutions. For example, the formulas (12. 28 p . (I2.29|) and (I2.30P 
correspond to (26), (25) and (30) of [3]. 

The regularity of the solutions and/or the positivity of the resulting weight functions 
etc. must be verified in each specific case for the chosen parameter values. This is in good 
contrast with the generic multi- indexed Laguerre and Jacobi polynomials in [3], in which the 
conditions for the parameters are given (eqs.(23) and (24) of [3]). 

When there is no singularity in rj G (— 1, 1), the norm of the jV-th transformed solution 
<f>^ of the n-th the eigenfunction (j) n {x) is related to that of the original <f) n (x): 

M N L 

j=i j=\ j=i 

where £j is the energy of the j-th seed solution. 

3 Higher degree of apparent singularities 

We are going to investigate the relationship between the zero of the Wronskian and the 
singularity of the deformed Schrodinger equation. 
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Let us start with a general Schrodinger equation (12. ip with an unspecified potential 
U(x). Suppose Xq, a regular point of the original potential U(x), is a zero of the Wronskian 
W[(fi, . . . , (pfj}(x) (12. 4p with multiplicity m, i.e. 

W[^, . . . , <^](x) = c (x - x ) m + 0(x - x ) m+1 , m G Z >0 . (3.1) 

Then it is a singularity of the potential of the deformed Hamiltonian such that 

x « a; , U^(x) = U(x ) + - — ^ m + regular terms. (3.2) 

{X - Xq) Z 

The point x is an apparent singularity of the deformed Schrodinger equation, because the 
monodromy is preserved by Darboux transformations. The characteristic exponents at the 
apparent singularity x Q are determined by 

p(p - 1) - 2m = 0, (3.3) 

and the difference of the exponents is an integer because of the apparency of the singularity. 
Combining it with that m G Z>o, the admitted values are 

m = l(l + l)/2, p = -l, l + l, (ZeZ >0 ). (3.4) 

Hence the degrees of zeros of the Wronskian are restricted to m — 1, 3, 6, 10, ... . In the case 
m = 1, the exponents are —1 and 2, and the higher exponents of the singularity are realised 

only in the cases m — 3, 6, 10, 

Next we investigate the condition that the Schrodinger equation of the deformed Poschl- 
Teller potential (I2.27P has an extra singularity in t] with higher exponents on a curve g = g(h) 
(or equivalently, h = h(g)) in the two dimensional parameter plane (g, h). This means that 
the polynomial T>(rj; g(h),h) has a multiple zero except for 77 = ±1. Hence a necessary 
condition for the curve g = g(h) is that the discriminant of the polynomial T>(rj; g(h), h) in 
the variable rj vanishes. It is not sufficient because the multiple zeros of T>{r\\ g(h),h) may 
be realised at 77 = ±1 or the coefficient of the leading term of T>(rj; g(h), h) may be zero. If 
the polynomial T>(rj; g(h),h) has a multiple zero at 77 = i]o(h), then the degree of the zeros 
are restricted to m = 3, 6, 10, ... . By fine-tuning the parameters on the curve g = g(h), the 
confluence of some apparent singularities is achieved. The fine-tuning has another effect. At 
the multiple zeros, the numerator polynomial Q(rj) and the denominator polynomial T>(rj) 
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are no longer relatively prime. While the denominator T>(rj) has an m = 1(1 + l)/2 zero, the 
numerator polynomial Q(rj) has an 1(1 + l)/2 — I = 1(1 — l)/2 zero: 

V( V ; g(h), h) = ( V - m (h)) l ^ 2 w(m 9(h), h), (3.5) 
w(r];g(h),h) = (ri-ri (h)) 1 x q(r]), deg(g) —£ — l, (3.6) 
Q»(»7; = (V - Vo(h)) l{l - 1)/2 V n ( m g(h), h), V n ( Vo (h); g(h), h) ^ 0. (3.7) 

The function (ff\x) = (l-r]Y^I 2 (l+r]f^' 2 V n (T]]g(h),h)/w(T]]g(h),h) defined by f[2~29l 
satisfies the deformed Schrodinger equation and the exponents are —I and I + 1. If I = 2, 
the denominator Wronskian has a triple zero, and the degree one polynomial is cancelled by 
that of the numerator. It should be stressed that the normalisation of the polynomials {V n } 
is uniquely determined by that of 10(77; 9(h), h). 

Hereafter, for simplicity of presentation, we simply write V n [r\\ g, h), etc. when the curve 
g = g(h) (or h = h(g)) is pre-specified. 

4 Four groups of new Jacobi systems 

Now we present several concrete examples of interesting systems exhibiting apparent singu- 
larities with the characteristic exponent p = —2. We shall be interested in new quantum 
systems generated by the Darboux transformation using the following eight sets of seed 
functions: 

[g, 4>™}; [ft, 4>™}; [ft, ft, ft 1 ]; [ft, ft\ ft 1 ]. 

The four sets in the second row are obtained by replacing in the first row the type I seed 
functions by the type II ones. Systems generated by other choices of the set of seed functions 
can be studied accordingly by the method described in the main text. 

4.1 The four groups of Jacobi systems 

We look for the situation where the polynomial part of the Wronskian of the set of seed 
functions, T>(rf) (I2.23p . contains a cubic zero outside of the physical domain x G (0,7r/2) 
(or 77 G (—1,1)). It turns out that such a situation exists for these eight cases above when 
g = g(h) is a linear fractional function. Note that it is equally well to consider h = h(g) as a 
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function of g in each case. But taking g as a function of h allows us to consider the h — > ±00 
limit, which will take the Poschl- Teller potential to the radial oscillator considered in jl]. 



It is discovered out that these eight cases fall into four groups of two members each, 
such that one member of a group can be obtained from the other by a shift in the value of 
h, which is determined by the universal formula for g and h, (I2.28P derived in section 12.31 
This is because, the main part of the denominator polynomial T>(rj] g(h), h), which has a 
cubic zero, is the same for the two members when expressed in terms of h. Moreover, the 
polynomial up to a multiplicative constant is identical for the Groups I and II and it has 
the same form for the Groups III and IV, with only the value of a parameter 7 is different. 
This means that the deformed potential U^(x) is identical for the Groups I and II, and it 
has the same form in the Groups III and IV, when expressed in terms of g and h. 



• Group I: (a) g = 5g=|l, g = g, h = h - 2; 

Group I: (b) [%,%,jF\,g=*g^, g = g-3, h = h + l. 
For both cases, g = ^rEi , g + h + 2 = g + h, (L — 1). 

• Group II: (a) [${, 0™], g = g = g, h = h - 2; 
Group II: (b) 9 = ^ 9 = 9~2, h = h. 
For both cases, g = g + h + 2 = g + h, (L = 1). 

For Groups I & II: T>i(r]; g, h) = (Zhrj - 2<q - 2h - i f x constant, (4.1) 




2hr] -277-2/1-I (2h7] - 2r] - 2h - l) 2 



(4.2) 



• Group III: (a) [^ 2 I ,0? 1 ], 
Group III: (b) [4>\, fe, 4>? 1 } 
For both cases, g = 




9 = ^^-1 9 = 9-2, h = h; 
9 = 9 = 9 + 1, h = h-3. 

g + h + 2 = g + h, (L = 1). 



• Group IV: (a) [<, <$ n ] 
Group IV: (b) [^i,^ 11 ,^ 11 ]. 



For both cases, g = 




For Groups III & IV: V m (r); g, h) = (2/177 - 2/i + 3) 3 x constant, 



(4.3) 
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< , w = M^ i ) + 2^ i ) _ + - + 
I — 77 1 + 77 

3(2/t-3) 9(4/i-3) 



(4 4) 

2^77 — 2/i + 3 (2/ir/-2/i + 3) 2 _ ' 
7 = 2 (for Group III), 4 (for Group IV). (4.5) 

The normalisation of the (fine-tuned) multi- indexed polynomials {V n (r])} is fixed by 

<p^\x) = (1 - r,yw\i + v ) m/2 r n (m g(h), h)/w( m g (h), h), (4.6) 

w(f]) = Wi(i]; g, h) = f (2hr] — 2r] — 2h — l) 2 , for Groups I & II, 
w (rj) = w 2 (m 9, h) = (2hr] -2h + 3) 2 , for Groups III & IV. 

For these to be non-singular in the physical domain 77 G (—1,1), one must restrict 

h > - and h ^ 1, : I, II, h< ^ and h ^ : III, IV. (4.7) 

The conditions h ^ 1,0 are also needed, otherwise the denominator in each case will be a 
constant, which is not the case we want to consider here. The square integrability of the 
eigenfunctions requires 

9>~\, ~ h >~\- (4-8) 

In Appendix we list the possible range of h for each specific case. 

We note here that g\ = 1 — gn and g~m = l — giy. This explains why the potentials in r\ are 
the same for Groups I and II, and for III and IV. By equating h for the two cases in the same 
group, one obtains a transformation between the two systems. Hence by setting h — > h + 3 in 
the potential, the eigenvalues, the eigenfunctions etc in Case la (with g expressed in terms 
of h), one obtains the corresponding formulas for Case lb. Similarly, setting h — > h + 2, 
h — > h — 3 and h — > h — 2 map Cases Ha to lib, Ilia to Illb and IVa to IVb, respectively. 
These are the required shiftings mentioned above. 

Here are the Fuchsian differential equations for the rational functions y n (r]) = f V n (77) jw {jf) , 
n = 0, 1, ... , with a p = —2 apparent singularity. They are obtained by rewriting the 
Schrodinger equation with the potential U^(x) f)4.2|) . (14.41) : 

(1 - v 2 )y'M + (h-g-(g + h + 1)77) y' n (v) + { + 1)0 + 9 + h + 1) 

+ 6 ^ + 1 ) + 18(4 "- 1} \v(v) = (49) 

+ 2^- 277 - 2^-1 + (2/177 -2 V -2h- 1)2 J MJ U ' 
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for Groups I and II, and 



(1 - v 2 )Vn(v) + (h-g-(g + h + l) v ) y' n { v ) + { (n + |)(n + g + h + |) 

6(2/i-3) 18(4/i-3) 1 




(4.10) 



for Groups III and IV. In these equations, g = g{h) should be replaced by its respective value 
for each group as given before. These are second order differential equations with regular 
singularities at rj = ±1, oo and rj = (2h + l)/2(h — 1) (for Groups I and II), or 1 — 3/2h (for 
Groups III, IV) with the conditions (14. 7p . 
The polynomials V n satisfy 



- {3(7 + 2n) 2 + 2h [3(7 + 2n) - I677 + (7 + 2n) 2 (r/ - 1)] + Ah 2 (A + 7 + 2n)(r) - 1) 

+2g [3(7 + 2n) + 2h{A{ V + 1) + (7 + 2n)( V - 1))] } V n ( V ) 
+4 {(2h 2 (7] - l) 2 + 3r] + g [3 + 2h{r) - 1)] ( V + 1) - 7i(6?7 + 5) (77 - l)}V n (v) 



for Groups III and IV. 

Appendix |A] presents the essential data for these new polynomials for one case from each 
group. The other member can be obtained by the shifting mentioned in the last paragraph. 
Examples lib and Ilia were discussed in [14J from a different point of view. 

4.2 Heun's systems 

All the examples in the previous subsection have four regular singularities, three at finite 
locations and one at 00. Thus they belong to the Heun's systems with one apparent singu- 
larity of characteristic exponent p = —2. Let us discuss the Heun's properties by taking the 



{l + 2n + n 2 - 2/i 2 (3 + n)(r] - 1) - 677 + 4m? + 2n 2 rj 



+ h[-l + n(5 - 2r]) - 2n 2 (r] - 1) + 12r/] + g [5 + n + 6ry + 2nr] 
-2h(l + n(r] - 1) + 377)] } V n {rj) + {-8 + 2h 2 ( V - l) 2 - 77 + 6r/ 2 
+g(-l + 2h{ J] - 1) - 2/7)(r7 + 1) - ^(77 - 1) (877 + 9)} ^(77) 
+ (-1 + 2/1(77 " 1) " 2r/)(r/ 2 - V)V'M = 0, 



(4.11) 



for Groups I and II, and 



+ [3 + 2/1(77-1)] (77 s - 1)^(77) = 



(4.12) 
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example, Group I (a) with g = 3(h — 3) /{Ah — 9). The T>(r]) polynomial (14. ip in terms of 
h — h + 2 reads 

T>(r]; g, h) = (2(h - 3)77 - 2h + 3) 3 x constant. 

By the change of the independent variable 77 = 1 — 2z, the equation for the two-indexed 
Jacobi polynomial y n = V n (r]; 3(h — 3)/(4h — 9), h) is rewritten into the Heun's differential 
equation (11. ip : 

d 2 y n , (g + l/2 , h - 3/2 I \ dy„ 



r, qU. 1 n 1 f, QU 3(4fe-9)» 2 +3(4fe 2 -6/i-9)n+2fe(10fe-27) 

-n 6)[n + g + n 6)z 4(h-3)(4h-9) 



y„ = 0, n = 0, 1, 



(4.13) 

The exponents at z = 0, 1 and the position of the extra singularity are specialized by the 
parameter h, and the exponents at z = 00 (i.e. — n — 3 and n + g + h — 3) depend on the 
discrete parameter n. The accessory parameter (q in equation (II. ip ) is specialized so that 
the extra singularity is apparent and the equation has a polynomial solution. As shown in 
IA.lt the equation has an extra solution for n = —2, (I A . 2 [) . 

It is interesting to note that the accessory parameter q in (14.131) has a very special 
dependence on a, (3 (the exponents at 00) and the location t of the extra singularity: 
3(4/i - 9)n 2 + 3(4/i 2 - 6h - 9)n + 2h(10h - 27) de f 



4(/i-3)(4/i-9) 



q = a(3t + q r 



def ( , o\ /P def , 1 t o , def 3 dcf 2h - 9 

a = — (ra + 3), p = ri + (7 + /i — d, r = —-77 -r, qv 



4(/i-3)' ^ 2(/i-3)' 
Here, the reduced accessory parameter q r does not depend on n. This simplifies (I4.13P to 
d 2 y n fg + 1/2 /i - 3/2 4 \ 4 , / a/3 g r 



dz 2 + V z + z-1 z-t) dz + \z(z-l) z{z- l)(z-t) ' ' /,! °' 



For all the examples in section 14.11 we have 

q = a(3t + q r , q r : does not depend on n, (4-14) 
and the following Proposition. 

Proposition 4.1 The polynomial solutions of the eight examples in section I4TTJ can be writ- 
ten in the following special form of Heun's equation with n-independent q r : 

d ^ + (i + A 4 -) d ^+( «p — * — \ =0j (415) 



dz 2 \z z — 1 z — t)dz \z(z — 1) z(z — l)(z — t) 

13 



q r = 



(2 + g + h-4gh)/{2_(h-l)} Groups I & II, 
(-6 + g(3 - Ah) + 7h) /2h Groups III & IV, 



1 + S = a + f3 + 5. (4.16) 



This can be verified in each case. The above relation (I4.14p can be rephrased that the pole 
residues of the extra singularities are n-independent. This is also the case for the Xg Jacobi 
polynomials which have p = — 1. See (3.1 1)-(3.12) in [15]. The same statement applies to 
the generic multi-indexed Jacobi polynomials, (4.6)-(4.7) in [3]. 

Another example is the case [</>!, with g = 1 — h. The denominator polynomial is 
simply 



The degree of the two- indexed Jacobi polynomial V n {r\\ 1 — h,h) is n + 2, and the function 
(1 — ry)( 1_/l )/ 2 (l + r]) h / 2 Vn(Vi 1 — h, h)/r] 2 is a solution of the deformed Schrodinger equation. 
In this case, however, the extra singularity at rj = is located in the physical domain, 
and it destroys the self-adjointness of the Hamiltonian as well as the quantum mechanical 
interpretation. From a viewpoint of mathematics, this case is out of the Sturm-Liouville 
theory, but it satisfies a second-order differential equation. Namely, the two-indexed Jacobi 
polynomial y n = V n {r\\ 1 — h, h) satisfies 



and it also satisfies Heun's differential equation by setting r\ — 1 — 2z. 

Similar phenomena happen in the cases = 1 — h}, {[0?, ^3], g = 1 — h}, 



Note that the Darboux-Crum transformation was already applied to Heun's equation in 
|16j . However the motivation and the results in the present paper are different from those 
in PS]. 

4.3 Limits to Laguerre systems 

We now discuss the relations of these Jacobi systems with the Laguerre systems presented 
in [I]. It is well-known that by a simple limiting procedure, one can obtain the (deformed)- 
radial oscillator potential from the (deformed) trigonometric PT potential [17]. For this one 
makes use of the limit formulas of the base polynomials 



h, h) = rf x constant. 




/i + 3/2 /i-3/2 4\ dy, 
i] — 1 i] + 1 rj)dr l 




{[<f>? I ,<t> I i I ),9 = h}. 



[3— too 



lim P^' ±/3) (l - 2xf3~ l ) = L^\±x). 



(4.17) 
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By taking the limit of infinite coupling h — > oo with the rescaling of the coordinate: 

x = —=, < x < - ^> < x L < -Vh, (4.18) 
v/i 2 2 

we have, with i](x) = cos(2x) and r] L (x L ) = x l2 , the relation 

„( x ) = l-^ + 0(/T 2 ). (4.19) 

Taking the limit h — > ±oo (or equivalently, h — > ±oo), where the upper (lower) sign 
applies to Groups I and II (III and IV), we have 

lim h~ 1 £? l (g,h)=£^g) = 4n, (4.20) 

h— »±oo 



and 



lim 

h— »±oo 



l 3 48 288 / 1'. 

* -I? + 4^ T (i?^- 2 r 7+ 2 ' (l ' 21 



with 7 = 2 for Groups I to III, and 7 = 4 for Group IV. These are the potentials of the 
Laguerre systems presented in |4], where g = g(h) takes on the respective value for each 
group at the limits h — > ±oo, namely 

g^\ (II, III), ~ (I, IV). (4.22) 

The Group I corresponds to the case (A), II to (B) and III to (D) in [4]. 
The equations (14.91) and (14.101) for 7/(77) reduce to 

vv" + (9 + \- v)y' +( n + l~ ^±1^ ) y = °' n = > 1 >-> ( 4 - 23 ) 

for Groups I and II (upper sign and 7 = 2) and groups III and IV (lower sign and 7 = 2 and 
4, respectively). This corresponds to the equation (3.40), (3.49) and (3.55) in |4]. 

In the same limit, the polynomial V n reduce to the corresponding polynomial C n , which 
satisfies, from the limits of f)4.1ip and (14.121) for y(i]): 

477(477 + 3)^(77) - 2 (-3 - 6# + 3477 - 8gr] + 8t7 2 ) £^(77) 

+ 4 (11 - 8g + 3n + 4(3 + 77)77) £ n {v) = 0, (4.24) 

for Groups I and II, and 

477(477 - 3)^(77) - 2 (3 - 2^(47/ - 3) + 22r/ + 8r/ 2 ) C'M 
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+ (32 - 32g - 67 - Yin + 8(4 + 7 + 2n)n) C n (rj) = 0, 



(4.25) 



for Groups III (7 = 2) and IV (7 = 4). These correspond to the equation (3.41), (3.50) and 
(3.56) in 0]. 

The above results indicate that the systems described in Group I, II and III ( the potential, 
eigenfunctions, eigenvalues, etc.) reduce to the Cases (A), (B) and (D) Laguerre systems 
given in [4], respectively (for details of the Laguerre systems, see jl]). Group IV corresponds 
to a new Laguerre system. Since this is not discussed in [1], we shall give some details of 
this case in the Appendix. 

Specifically, the forms of U^\x) in (14. 2 p and (14. 4p imply that the eigenfunctions have 
the form: 

^x;S,/0 = - \s " V n (mg,h), 2 = 1,2. (4.26) 

Under the above limiting process, the prefactor of V n (r];g,h) (times the factor h~z) tends 
to 

T 3. —B— 

1 (4.27) 



(4r/ L ±3) 2 ' 

with g — 3/4 for Groups I and IV, and g — 1/4 for Groups II and III. The polynomial part 
h~ m V n (i]; g, h) (m = 2 for Groups I, II, and III, and m = 3 for Group IV) approaches to 

£ A,B,D,New ( ^ )? (42g) 

with g = 3/4, 1/4, 9/4 and 13/4, respectively. 

5 Other examples of higher degrees of apparent singu- 
larities 

5.1 The case [02 , 02 ] an d related cases 

In the case [(p 1 ^,^ 1 }, the Wronskian is written as Wf^ 1 ) 4 ) 2 1 ) { x ) — (1 — ri) 3 / 2 ~ g T)(ri), where 
T>(rj) is a polynomial of degree 4. The discriminant of the polynomial T>(rf) is written as 

-256(2/i + l)(2h - 3)(2h + 3) 2 (2h - 5) 2 (2g - 3)(2g - ll) 3 (2g - 5) 4 
x(g + h-4){g-h- 3){2gh 2 - 2gh - 3g - 5h 2 + 5h + 12) 2 . 
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We set g = (5/i 2 -5/i-12)/(2/i 2 -2/i-3), which is a root oi2gh 2 -2gh-3g-5h 2 + 5h+12 = 0. 
Then the polynomial T>(t]) is written as 

{(h + 2){h - 3)77 - {h 2 - h + 3)}{h(h - l)rj -{h 2 -h- 3)} 3 , (5.1) 

and the potential is 

2(g -2){g-3) + 2h(h - 1) _ + 2 + 24/i(/i - l){(/i 2 - - % - fc(/i - 1)} 



r] r] + l ^ ' (h(h - 1)77 - (h 2 - h - 3)) 



2 



8 (fr + 2)(fe - 3){(fr 2 - h + 3) V - (ft + 2){h - 3) 

((/i + 2)(/i-3)77- (h 2 -h + 3)) 2 ' l ' J 

Note that 77 = 1 — 3/{h(h — 1)} is an apparent singularity of exponent —2 and 77 = 1 + 

9/{(/i + 2)(/i - 3)} is that of exponent -1. 

The four cases <g n ], & 02 + [02, [0" 0" 2 U ] can be treated simi- 

larly to Groups I and II in section H] and the four cases r - ! 1 11 " " l| 1 11 111 1)1 



J 2iV2 h IVl ) 72 ) V2 J > LV2 5 VI )V2 J, 



Lvij V2j vi '^2 J corresponds to Groups III and IV. 
5.2 The case [<p\, ^i 1 ] and related cases 

In the case [<f)\, the Wronskian W[0J,\ (f) l 2 u }(x) = V(rj) is a polynomial of degree 4. The 
discriminant of the polynomial T>(rj) is written as 

-1024(2/i + l)(2h - 3)(2h + 3) 2 (2# + l)(2g - 3)(2g - 5) 2 (g - h - 3)(g - h - l) 3 

x (4g 2 h - Agh 2 - 30# 2 - 56gh - 26/i 2 + 69g + A3h - 12) 2 . 

The equation Ag 2 h - Agh 2 - 30g 2 - 56gh - 26h 2 + Q9g + A3h - 12 = is quadratic in both 
g and h. We find some lower lying rational solutions, (g,h) = (131/2,25/2), (115/2,16), 
(113/2,41/2). For these points in the parameter plane, the polynomial T>(r]) is proportional 
to 

(13 + 7r?) (5r? + 7) 3 , (45r/ + 107)(3t7 + 5) 3 , (77 + 3) (77 + 2) 3 , (5.3) 
respectively. The potential for the (g,h) = (113/2,41/2) case is 

def l^lll 41-39 /113 41\ 2 24(2r/ + l) 8(377 + 1) 
KU 2(1-77) 2(1 + 77) V 2 2 (t/ + 2) 2 ^ (t/ + 3) 2 • 

Note that 77 = —2 is an apparent singularity of exponent —2, 77 = —3 is that of exponent —1, 
and they are outside of [—1,1]. The corresponding (fine-tuned) multi-indexed polynomials 
{Pn(v)}i satisfy 

- (1 - V 2 )( V + 2) (17 + 3)^4^ + 2(116 + 32777 + 205^ + 3Q V 3 )^^- 

ar] z drj 
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+ (-2(146 + 355?/ + lllr/ 2 ) - An{n + 77) (77 + 2)(t? + 3)) P n (r?) = 0, (5.5) 
and the solutions are 

^(77) = -82 { [9(1517 + 1952r? + 825r/ 2 + lUr] 3 ) + n(n + 77) (2 + r/) 2 (3 + r/)] P^ w) (r]) 

+ (1 - ?? 2 )(2 + 77)(11 + 4r / )9 r ,Pi 56 ' 20 )(r / )} . (5.6) 

They satisfy the orthonormality relations 

- 1 (l-7/) 56 (l+7/) 20 



l + 2)4(^ + 3)^ 



V n (ri)Vm(ri)dri = H n 5 nm , (5.7) 



def 4 Mrc + 77) + 1102)(n(w + 77) + 1242)(n + l) 20 (n + 57) 20 

(2n + 77) ' 1 ' J 

The other parameter choices (g,h) = (131/2,25/2), (115/2,16) show equally interesting 
characteristics. The cases [03,0™], [</>!, 03 1 ], w£, <fI\ Wii ^ 1 ] have similar features. 

6 Summary and Comments 

Several families of infinitely many rational solutions are constructed of certain second order 
Fuchsian differential equations having four or more regular singularities. All the extra sin- 
gularities are apparent and their characteristic exponents are —2 or —1. They are obtained 
by fine-tuning the parameters of multi-indexed Jacobi polynomials [3] . In other words, these 
Fuchsian equations are derived from the Schrodinger equation with the Poschl- Teller poten- 
tial through Darboux-Crum transformations with the seed solutions having a cubic zero. The 
fine-tuning of the parameters (g, h), is realised as a curve g = g(h) (or h = h(g)) obtained 
from a particular root of the discriminant of the seed solutions corresponding to a cubic zero. 
These families of infinite solutions constitute orthogonal polynomials over (—1, 1) with the 
weight function of the form (1 — x) a (l + x)" / {(ax + b) A q(x) 2 } , in which q(x) is a polynomial 
in x. Being the special cases of the multi-indexed Jacobi polynomials, these polynomials 
have "gaps" in the degrees and they do not satisfy the three term recurrence relations. The 
eight examples presented in section four are divided into four groups, each consisting of two 
members. Reflecting the structure of the theory, the orthogonal polynomials of one member 
of the group, can be obtained from those of the other's by a simple shift of the parameters. 
Examples lib and Ilia were discussed in [14J from a different point of view. Properties of 
Heun's equations with an apparent singularity are discussed in some detail for two explicit 
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examples. By the standard limiting procedures leading from the Jacobi to the Laguerre 
polynomials, these special cases of multi-indexed Jacobi polynomials are shown to reduce to 
the special cases of multi-indexed Laguerre polynomials. Three of them had been explicitly 
constructed in [I]. 

It is a good challenge to construct solutions with higher degrees (p = —3, —4) of apparent 
singularities within the multi-indexed Jacobi polynomials. In the preparation of this article, 
we tried in vain some combinations of lower degree virtual state solutions. With increasing 
complexity in dealing with higher degrees of virtual state solutions, some new ideas might 
be needed. 
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Appendix 



A Basic data for four cases in section [4 
A.l Case la 

• 9 = 9 = 9, h = h - 2. 

• Explicit form of wave functions: 

V l *{rii h) = ( 1053 - 612/i + 96h 2 - 248/i 3 + 176/i 4 - 32h 5 - 162n + 108/m + lUh 2 n 
- lUh 3 n + 32h A n - 162n 2 + 288/m 2 - 168/iV + 32h 3 n 2 
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+ (1944 - 1026/i - 1740/i 2 + 1976/7 3 - 736/i 4 + 9Qh 5 + 48677 - 108/m 

- 432/i 2 n + 240/i 3 n - 32h 4 n + 486n 2 - 648hn 2 + 264h 2 n 2 - 32h 3 n 2 )r] 
+ (1458 - 4914/i + 6132/i 2 - 3520/7 3 + 944/i 4 - 96h 5 - 216/m 

- 12h 2 n + 144/7 3 77 - 32/7 4 77 - 216/m 2 + 168/iV - 32h 3 n 2 )r] 2 

+ (-1620 + 4248/i - 4020/i 2 + 1792/7 3 - 384/i 4 + 32/i 5 - 648n + 504h 2 n 

-2A0h 3 n + 32h A n - 6A8n 2 + 720hn 2 - 264h 2 n 2 + 32/i 3 n 2 )r/ 3 ) P^' 1 ' 2 ^- 1 ^^) 
+ (-9 + 4/i)(-l + ?7 2 )(3 - 2/i - 677 + 2/177) 

x (27 - 4/i 2 + 54r/ - 30/ir; + 4h 2 ri)d v Pi 9 - 1/2 > h - 1/2) (ri). (A.l) 

From (14.71) . one has h > 9/4 and h ^ 3. The square integrability (14. 8 j) requires a stronger 
condition h > 27/10. 

• Additional mode corresponds to n = —2 with 

P^n; h) = 27 - 4/i 2 + 2(/i - 3) (2/i - 9)77, (A.2) 

= _ g(9 _ uh + Ah 2 )/{Ah _ 9)= £lH ( gz| ) (A . 3) 

• Equation satisfied by V£'- 

(-9 + 4/i)(3 -2h + 2(-3 + %)(1 - t/ 2 )?/"^) 

- (27 + 6O/7 - 58/i 2 + 12/7 3 - 2(— 81 + 102/i - 48/i 2 + 8/i 3 )r/ 

+2(-54 + 57h - 19h 2 + 2/i 3 )(2r/ 2 - 1)) 1/(77) 

+ ( - 8/i 3 (3 + n) - 8/i 2 (-20 - 3n + n 2 ) - 27(-12 + 77 + n 2 ) 

+6/i(-63 + 5t7 2 ) + 2(-3 + h)(3 + n)(18 - I8/1 + 4/7 2 - 9ti + 4/777)77) 2/(^7) = 0. (A.4) 

• Orthonormality relations: 

r 1 n _ tiIs-V^i a. ft— 5/2 
jf^ \ 2/ ^V+3)' ^ " } ^ h) = n, m = -2, 0, 1, 2, . . . , 

^ = 2^ /l - 1 (4/7 - 9) 2 (f n - £t) (S n ~ Si 11 ) h n , n £ -2. (A.5) 

A.2 Case Ha 

•9 = Wig, 9 = 9,h = h- 2. 

• Explicit form of wave functions: 

, Ila, h) = (1 ~ + V)^ 1 IIa( . n = ()12 

^ n ^' j (2/777 -677- 2/7 + 3) 2 ' n [V ' 9h 

20 



V^im h) = {1458 - 918/i - 144/i 3 + I6O/1 4 - 32/i 5 - 144/m + 264h 2 n - 160h 3 n 
+ 32h 4 n - 162n 2 + 288/m 2 - 168/iV + 32/i 3 n 2 + (2187 - 2106/i 

- 936/i 2 + 1776/i 3 - 720/i 4 + 96/i 5 + 432/m - 600/i 2 n + 256/i 3 n 

- 32/i 4 n + 486n 2 - 648/m 2 + 264/iV - 32h 3 n 2 )r] 

+ (-3888/i + 6120/i 2 - 3624/r 5 + 960/i 4 - 96h 5 - 192h 2 n + 160/i 3 n 

- 32h 4 n - 216hn 2 + 168h 2 n 2 - 32h 3 n 2 )r] 2 + (-2916 + 6048/t 

- 4932/i 2 + 1992/i 3 - 400/i 4 + 32/i 5 - 576/m + 672/i 2 n - 256h 3 n 
+32h A n - 648n 2 + 720/m 2 - 264/iV + 32/i 3 n 2 )r? 3 } P^~ 1/2 ' h ~ 1/2) (r]) 
+ (-9 + 4/i)(-l + r] 2 )(3 -2h-6r] + 2hrj) 

x (27 - Ah 2 + 54r/ - 30/u? + Ah 2 r])d v P^- 1/2 ' h ~ 1/2 \r]). (A.6) 

From ( 14. 7p . one has h> 9/4 and h ^ 3. The square integrability ( 14. 8 p is satisfied. 

• An additional mode corresponds to n — —3 with 

pn>_ 1 cn a _ 12(27 - 20h + 4/^) ~ nI fc 

• Equation satisfied by V^ivi h): 

(-9 + Ah)(3 -2h + 2(-3 + /i)?7)(l - 77V' M 
- ((135 - 198/i + 96h 2 - 16/i 3 )r? 

+ (-3 + h){-9 - 20h + 12h 2 + (45 - 28h + 4/i 2 )(2r7 2 - l)))y'(ry) 

+ (3 + n) (27 - 20h + 4/i 2 - 9n + 4/m) (3 - 2h - 677 + 2hrj)y{r]) = 0. (A.8) 

• Orthonormality relations: 
1 (\ _ ri) 9 ~ l / 2 (\ + n) h ~ 5/2 

1 **! (2^7? - 677 - 2^ + 3)4 P " a ^ ; ^ fc ) = n ' m = ~ 3 ' °' 1; 2 ' ' ' ' 

C = 2^(4/1 - 9) 2 (E n - £{) (S n - 4 m ) /i n , n ^ -3. (A.9) 

A. 3 Case Ilia 

• 9 = = 9-2,h = h. 

• Explicit form of wave functions: 

€ la (v,h) = { \2hrj-2h + $" ^(mh), n = 0,1,2,...,. (A.10) 
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V^{ri; h) = {405 - 648/i + 72h 2 + 72h 3 + 80/i 4 + 32h 5 - 486n + 1620/m - 1728h 2 n 
+ 432h 3 n + 288h A n - 128h 5 n - \Q2n 2 + 648/m 2 - 936/iV + 576/iV 

- 128/i 4 n 2 + (810 - 1566/i + 612/i 2 + 408/i 3 - 192/i 4 - 96/i 5 + 486n - 2268/m 
+ 3456/i 2 n - 1584/i 3 n - 480/i 4 n + 384/i 5 n + 162n 2 - 864/m 2 + 1656/iV 

- 1344/iV + 384/i 4 n 2 )r/ + (270/i + 36/i 2 - 528/i 3 + 144/i 4 + 96/i 5 + 648/m 

- 19Uh 2 n + 1584/i 3 n + 96/i 4 n - 384/i 5 n + 216/m 2 - 792/i 2 n 2 + 960/i 3 n 2 

- 384/i 4 n 2 )r? 2 + (36/i 2 + 48h 3 - 32/i 4 - 32h 5 + 216h 2 n - 432h 3 n + 96/i 4 n 
+128/i 5 n + 72h 2 n 2 - 192h 3 n 2 + 128h 4 n 2 ) V 3 } pfo-V2,A-i/2) ^ 

- 9(-3 + 4/i)(-l + r/ 2 )(3 - 2/i + 2hrf){3 - Ah 2 - 2h V + Ah 2 r])d v P^- 1/2 ' h - 1/2) (rj). 

(A.ll) 

From ( 14. 7ft . one has h < 3/4 and h ^ 0. The square integrability ( 14. 8ft requires —1/2 < /t. 

• Additional mode corresponds to n = —2 with 

P n 2 Ia (?7; h) = 3 - Ah 2 + 2h(-l + 2h)r], (A. 12) 

£ nia = _ 8( _ 3 _ 2 ^ + 4^2 )/( _ 3 + Ah) = ^(gzil, /,). (A .13) 

• Equation satisfied by V™ 3 "- 

(-3 + 4/i)(3 - 2/i + 2/177) (1 - r/ 2 y'(r/) 

+ {9 + 2h 2 - 12h 3 + 2(9 - 6h - 8h 2 + 8h 3 )r] - 2h{3 -7h + 2h 2 ){2rj 2 - l)}y'{v) 
+ {-8h 3 (3 + n)- 8/i 2 (-3 + n 2 ) - 9(2 + 3n + n 2 ) + 6/i(l + 6n + 3n 2 ) 

+2/i(3 + n)(4/i 2 - 3n + /i(-6 + 4n))r/} y(rj) = 0. (A.14) 

• Orthonormality relations: 

^Illa def 29+ h-l {Ah _ 3) 4 ^ _ gll^j ^ _ glll^ n ^ _ 2 (A15) 



A. 4 Case IVa 

•9 = 3( 4fc-T 7 8) » g = g-3,h = h-l. 
• Explicit form of wave functions: 



, iv af ^ = (l-7y)V(l + iy)V 

^ n ^' j (2/ir/-2r/-2/i + 5) 2 7 n [V,9h ' ' ' 
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V™ a (r); ft) = { 133407 - 291870ft + 189396ft 2 + 30456ft 3 - 98352ft 4 + 51552ft 5 - 12096ft 6 
+ 185640™ - 460616/m + 440564ft 2 n - 191368ft 3 n + 24144ft 4 n + 10336ft 5 n 

- 4352ft 6 n + 54145n 2 - 147238ftn 2 + 165212ftV - 98712ft 3 n 2 + 33568ftV 

- 6272ftV + 512ftV + 128ft 7 (9 + 4n) + [-194400 + 574560ft - 581760ft 2 
+ 141696ft 3 + 160128ft 4 - 137088ft 5 + 41472ft 6 - 164976n + 604208/m 

- 830840ft 2 n + 526768ft 3 n - 129120ft 4 n - 16384ft 5 n + 14336ft 6 n - 48118n 2 
+ 187684/m 2 - 288488ftV + 227232ft 3 ™ 2 - 97792ft 4 n 2 + 22016ft 5 n 2 

- 2048ft 6 n 2 - 512ft 7 (9 + 4n)]r? + [-97686 + 4860ft + 300672ft 2 - 256608ft 3 

- 38016ft 4 + 124416ft 5 - 51840ft 6 - 67032n - 77448/m + 449940ft 2 n 

- 490632ft 3 ™ + 194688ft 4 - 2496ft 5 n - 16896ft 6 n - 19551n 2 - 17934/m 2 
+ 134904ft 2 n 2 - 176112ft 3 n 2 + 102336ft 4 n 2 - 28416ft 5 n 2 + 3072ft 6 n 2 

+ 768ft 7 (9 + An)]rf + [32400 - 23760ft - 102240ft 2 + 147744ft 3 - 33408ft 4 

- 43776ft 5 + 27648ft 6 + 49728™ - 71872/m - 67232ft 2 n + 179392ft 3 n 

- 108960ft 4 n + 12800ft 5 n + 8192ft 6 n + 14504n 2 - 24416ftn 2 - 13352ft V 
+ 55008ft 3 n 2 - 45568ft 4 n 2 + 15872ftV - 2048ft 6 n 2 - 512ft 7 (9 + 4n)]r/ 3 
+ [-3240 + 4320ft + 9576ft 2 - 21168ft 3 + 9648ft 4 + 4896ft 5 - 5184ft 6 

- 3360n + 5728/m + 7568ft 2 n - 24160ft 3 n + 19248ft 4 n - 4256ft 5 n 

- 1280ft 6 n - 980n 2 + 1904/m 2 + 1724ftV - 7416ft 3 n 2 + 7456ftV 
-3200ft V + 512ft 6 n 2 + 128ft 7 (9 + 4n)]r? 4 } pjg-i/**-m {jfi 

+ (-7 + 4ft)(-l + rj 1 ) {336n(5 - 2r]f(r] - 1) + 98n 2 (5 - 2r]f(r] - 1) 

+ 16ft 4 [120 + n(38 - Urj) + 8n 2 (r] - 1) - 69r]](r] - l) 2 

+ 32ft 5 (9 + An)(r] - l) 3 + 3(-2048 + 4769r? - 2056r? 2 + 388r/ 3 ) - 8ft 3 (r/ - 1) 

x [-438 + 468?7 - 93r? 2 + 36n(2 - 7r] + Arf) + 8n 2 (17 - 28r/ + llr/ 2 )] 
+ 4ft 2 [9(-32 + 178?7 - 152ry 2 + 51r? 3 ) + 6n 2 (-143 + 333r? - 249ry 2 + 59r? 3 ) 
+ 2n(-953 + 2373r/ - 1896r/ 2 + 476r? 3 )] - 2ft[28n 2 (-85 + 174r? - lllr? 2 + 22r/ 3 ) 
+ 8n(-895 + 1863r/ - 1212r/ 2 + 244r? 3 ) 

+3(-1447 + 3859r/ - 2144r/ 2 + 488r? 3 )]} d v Pjf- 1/2 ' h - 1/2 \ri). (A. 16) 

From (14.71) . one has ft < 7/4 and ft ^ 1. The square integrability (14. 8 j) requires a stronger 
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condition 1/2 < h < 13/10. 

• There are two additional modes corresponding to n = —3 and n = —2: 

(i) n = -3: With 

P% = (2h-5)r,-2h + 2, g™» = - 12( ~ 3 + 2K) = h). (A.17) 

(ii) n = -2: With 

piva = 1Q7 _ g8/l + 24/l 2 _ 32/l 3 + 16/l 4 

- 8(-25 + 35/i - 14/i 3 + 4/i 4 )r/ + 4(-l + h) 2 (ll - I2h + 4/i 2 )r/ 2 ), (A.18) 

= - ^h-^ = ^Wfr* V (A.19) 

• Equation satisfied by P^ A (^; h): 

(-7 + 4/i)(5 -2h + 2(-l + %)(1 - 77 2 )2/"(^) 

+ {(-5 + 62/i - 64/i 2 + 16/i 3 )?7 - (-1 + h)(19 - 28h + 12h 2 + (5 - 2h) 2 {2rf - l))}y'{r]) 
+ {-(-5 + 2/i)[-12 - 24n - 7n 2 + 4/i 2 (4 + n) + 4/i(-4 + 2n + n 2 )] 

+2(-l + /i)(4 + 4/i 2 + 4/i(-2 + n) - 7n)(4 + n)r?} y(r/) = 0. (A.20) 

• Orthonormality relations: 

/ ? V (2hri -2 V -2h + 5y V "* im k) V ™ a{m k) = /l " Va<5nm ' n ' m = _3 ' ~ 2 ' °' 2 ' ' ' ' ' 
C = 2* +fc - s (4/i - 7) 6 (V n - £?) - A" 1 ) (£» - fcn, n ± -3, -2. (A.21) 

A. 4.1 Laguerre limit 

In the h — > — oo limit as described before, the polynomial P^ Va (?7; h) (n = 0, 1, 2, . . .) tends 
to 

£^ cw (r/) = {5265 + (-5616 - 2592n)r? + (11232 + 5184n)r/ 2 

+ (-6912 - 3072n)?7 3 + (2304 + 1024n)^ 4 } L^ 3/4) (rj) 
+ {1620?7 - 4(252 + UAn)rf - 4(-816 - 384n)^ 3 

-4(576 + 256n)^ 4 } d v L^ 3/4 \r]). (A.22) 

Here we have adopted rj instead of r] L for simplicity of presentation. The two additional 
modes approach to 

P w 3 a ->■ -(4r/ + 3), (A.23) 
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(A.24) 



They satisfy the equation 

4 V (4 V - 3)C*™(r))" - (Ar, + 3) (4^ + 5)£f w (r?)' 

+ 4(4m? + I677 - 3n - 4)£^ ew (r/) = 0, n = -3, -2, 0, 1, 2, . . . (A.25) 
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